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Abstract 



It is shown that the characteristics of the mesoscopic fluctuations in the con- 
f~^ ' ventional quantum-diffusion model and the model of the non-coherent ('clas- 



sical') diffusion in media with long-range correlated disorder are quite similar 
in the weak-disorder limit. The relative values of the variance and of the high- 
er ' order moments of the fluctuations in one model are obtained from those in 
another one by substituting a proper weak-disorder parameter. As behaviour 
of the ensemble-averaged diffusion coefficient is quite different in these models, 
it suggests that the mesoscopic properties, including a possible universality 



^ ' in the distribution functions, could be independent of the behaviour of the 



averaged quantities. 



I. INTRODUCTION 

The absence of self-averaging in mesoscopic conductors (see recent collections of review 
papers ||1]J^) is known to result from the interference effects in quantum disordered systems. 
The universal conductance fluctuations (UCF) in such systems are the most famous example 
of phenomena caused by the absence of the self-averaging. Conductance (or resistance) 
fluctuations are also characteristic of mesoscopic semiconductors with hopping conductivity. 
In this case, the absence of self-averaging has nothing to do with the quantum coherence, 
and is due to the fact that conductance is determined by the most favourable channel 
which should be totally different for different realizations of disorder, and could be easily 



changed by applying external fields. These fluctuations are by no means universal, and 
their amplitude could be much larger than an appropriate average value, in contrast to 
the UCF in the mesoscopic conductors whose amplitude is much smaller than the average 
conductance. In the absence of the quantum interference, transport properties are governed 
by the probabilities rather than by amplitudes of wave functions, so that the mesoscopic 
effects may be described in the framework of classical models. In general, mesoscopics of 
systems with and without quantum coherence are totally different. 

In this paper, however, we will consider some models of classical transport in disordered 
media 0- [§] whose "mesoscopic" properties show remarkable similarity [|^,|^ with those 
of the conductance fluctuations in metals ||TO|JTl|| . Firstly, the fractional value of the vari- 
ance of the fluctuations has the same dependence on the effective disorder parameter as 
that in the quantum diffusion problem, although the fluctuations in the former case are 
not universal, in contrast to the latter. Moreover, all the higher moments of the fluctua- 
tions have the same (nontrivial) dependence on the appropriate weak-disorder parameter. 
It makes the distribution functions of the fluctuating transport coefficients to have the same 
shape in the both cases: almost Gaussian in the bulk of the distribution but with slowly 
decreasing lognormal tails. Such a similarity is even more striking, taking into account that 
average values of the transport coefficients (conductance, diffusion constant, etc.) have abso- 
lutely different dependence on the disorder. The reason is that the "mesoscopic properties" 
could be unrelated to those particular features of a model which determine behaviour of 
the ensemble- averaged quantities {e.g. the conductance). Instead, they could be determined 
by some general symmetry that might appear independent of details of this or that model, 
and even of the presence or absence of the quantum coherence. We shall discuss this in the 
conclusion but now let us begin with describing the classical model under consideration. 



II. MODELS OF CLASSICAL DIFFUSION IN WEAKLY DISORDERED MEDIA 

First of all, to be able to compare transport properties in models of classical and quantum 
diffusion, we need to define a classical model for diffusion in a weakly disordered medium. In 
a generic quantum diffusion model (the Anderson model, or the model of free electrons in a 



random potential), the disorder parameter is directly related to the transport properties [jT2 
it is just the inverse dimensionless conductance, 1/g ~ l/{pplY~^, and the dimensionless 
conductance g is related to the diffusion coefficient D = VpTJd by the Einstein relation, 
g = Cdi'D {I is the mean free path, Cd is the constant depending only on the dimensionality 
d). In contrast to this, the parameter of quenched disorder in classical transport can be 
totally unrelated to the diffusion coefficient. We consider a generic model for non-coherent 
hopping transport, and discuss different ways of introducing a weak disorder into it. 



A. Lattice hopping models and continuum limit 

The lattice hopping model may be defined by the master equation for the probability 
Pr{t) of finding the particle at the r-th site at the moment t: 

-^ = Y1 {Wrr'Pr' " Wr'rPr) ■ (l) 

Here Wrr' is the probability of hopping from site r' to site r. We call this model classical just 
because transport is described in terms of the probabilities rather than amplitudes, although 
the origin of hopping could be quantum as well. 

In the continuous limit, this model goes over to the classical Fokker-Plank (FP) equation: 



d 



p(r,t) = 0, (2) 



where a is a vector index in the d-dimensional space. The parameters of this equation are 
related to the hopping probabilities W as follows: 

D{r) = ^-j:{r-rrWrr' (3) 

r' 

v(r)= Y.i^-r')Wrr' (4) 

r' 

On a regular lattice in the absence of external fields, the hopping probability W^ is sym- 
metric, W^r-r' = W^r'-r5 ^o that the drift velocity v, Eq. (§), vanishes and only the diffusion 
term survives in the FP equation (|^), with the diffusion coefficient D given by Eq. (^). 

Now one introduces a weak disorder into the model, making the hopping probability W 
slightly fiuctuating around its mean value W^ on the regular lattice, 

Wrr' = lVr_r' + ^Wrr' , (5) 

with 6Wrr' describing the quenched weak disorder on the lattice. It contributes to both 
D and v in the continuous limit of the model, Eq. (^. Were the hopping probabilities 
W = W'^ + 6W kept symmetric, v would vanish and the weak disorder would enter the 
model only via slightly fiuctuating diffusion coefficient D. We will show that such a weak 
disorder is irrelevant, in a sense that it does not effect large-scale (or long-time) properties 
of the hopping model. The only source of a relevant weak disorder could be in the drift 
velocity v, i.e. in the asymmetry of the hopping probabilities 6W. 

The random asymmetry that leads to the fluctuating v could be due to the presence 
of magnetic or charged impurities. Following ref. |]^, we will show that the presence of a 
random electric field E(r) of charged impurities leads to the potential random drifts (with 
curlv = 0) in the continuum-limit model, Eq. (0). If E is the only source of randomness, 
the probability of thermally activated hops at a distance b is 



Wr,r+b = W^exp 



eE(r)b' 
2kT 



(6) 



Assuming then a random distribution of the charged impurities on the lattice and a global 
electro-neutrality, and representing the Fourier transform of E(r) as 

E(q) = 2^ Cj exp(2q ■ r •), 

eq y 

with Cj = ±e being the charges of impurities with random coordinates r • and e being the 
dielectric constant, one finds 

{E^{q)Ep{q )) = 3" V (7) 

where C is the concentration of the charged impurities. Let us also assume for simplicity the 
nearest-neighbor hopping only, and take the weak-disorder limit which allows to expand the 
exponent in Eq. (|^). Then, substituting 6W oc E into Eq. (^) that define the drift velocity 
in the continuum limit, one arrives at the model described by the FP equation (0) with 
random drifts v with zero average and the correlation function 

Mr)vp{r')) = 7oFa/3(r - r') (8) 

where the Fourier transform of Fa/3{r — r') is given by ^^/^(q) = qaqp/q^- Obviously, 
curl V = so that we will refer to this model as one with the potential random drifts, 
or just the potential disorder model. The strength of disorder 70 in Eq.(P) is related to the 
parameters of the hopping model as 

^°~ £2(^2^)2 1^ 2 ) - e^ikry '' ^ ^ 

where z is the coordination number of the lattice, a is the lattice constant, Wq = Wj^ is 
a regular part of the nearest-neighbor hopping probability, Eq. (^, and Dq is the average 
diffusion coefficient. The randomness in W, both symmetric and asymmetric, results also 
in spatial variations of the diffusion coefficient, D{r) = Dq + SD{t). The variations can be 
described in terms of some distribution P{D) governed by the cumulants 

{{5D{v{)5D{v^))) =2!g(2)5(ri-r2), 

(10) 

((<5D(ri) . . . 5D{vs))) = s\ g(^)(5(ri -r^)... 5{t, - r,) , 

which may be directly expressed via the cumulants of 5W . The absence of spatial correlations 
of D{y) in the continuum limit results from taking into account next-neighbor hopping only. 
In the next section we show that it is the presence of the random drifts, Eq. (H), that 



changes drastically the long-time (or long-range) asymptotic properties of the transport 
coefficients while the disorder in the diffusion coefficients is totally irrelevant for the behavior 
of the averaged quantities. On the other hand, we will show that whatever were the initial 
distribution of the diffusion coefficient characterized by the couplings g^''^ Eq. (|TD|), its 
asymptotic characteristics would also be drastically changed due to the influence of the 
random drifts. 

Similarly, in the presence of a random magnetic field one arrives |^ at the model of 
solenoidal random drifts with zero average and the correlation function (^) with the Fourier 
transform of ^^^(r — r') given by FQ,^(q) = dap — qa'ipl'f'- Finally, one may consider a model 
with the short-range correlation function [Q, Fapij — r') = dapSir — r'), although the latter 
one is not directly related to the lattice hopping models (besides the case where a specific 
relation is imposed between an external magnetic field and a disorder parameter [J13l)- So we 



can consider the models described by the FP equation @) with the irrelevant weak disorder 
in the diffusion coefficient, Eq. (p!0|), and the relevant one in the drift term, Eq. (|]), which 
is characterized by three possible choices for the correlation function: 

g^g/j/g^, potential disorder (a) 

'ap = \ 5ap — QaQp/q'^, soleuoidal disorder (b) (H) 

6ap, mixed disorder (c) 

The FP description above is the basis for the further considerations. An appropriate 
Langevin description is also quite useful as it provides a clear qualitative picture. A diffusing 
particle can be considered as a random walker exposed to a strong thermal noise ('random 
winds') and a relatively weak quenched disorder ('random drifts'): 



V r 



r7(r,t) (12) 



In order to reproduce the FP equation (^, the thermal noise r} is chosen to be the Gaussian 
white noise with a zero average and 



r]a{r, t) r]p{r, f) = 2D{r) S^p 5(t - t'), (13) 

where D{r) is the weakly fluctuating diffusion coefficient with the distribution (plQl). Were 
the quenched disorder strong, the random walker would mainly follow the drift lines of the 
field V with rare hops from one line to another due to the presence of the random winds. 

We consider, however, only the case of the weak disorder when the motion is mainly 
governed by the random winds. In the absence of v, such a motion would be pure diffusive. 
The weak random drifts correct this diffusive motion as the random walker experiences also 
some almost ballistic motion along the drift lines, albeit short and frequently interrupted by 
hops between the lines induced by the random winds. The potential model is characterized 



by the presence of sinks and sources in the field fines. Witfiout tfie tfiermal noise, tfie 
particle would end at some sink (i.e. would be confined in a restricted region of space). 
Then it is natural to expect that the infiuence of the weak potential drifts lead to sub- 
diffusion behaviour in the long-time limit. In the solenoidal field, there are only closed drift 
lines whose presence could help the particle diffusion, so that one could expect a super- 
diffusion behaviour. And in the case of the mixed disorder, these two trends could partially 
cancel each other. We show further how to confirm this qualitative picture. 



B. The effective functional 

A standard way to describe the long-range and long-time asymptotic processes is to 
derive and solve the renormalization group (RG) equations for the quantities of interest. To 
this end, we introduce the effective field-theoretical functional describing the properties of 
the model. 

All the transport properties may be expressed in terms of the Green's function of the FP 
equation (^. In the absence of the drift fields, the Green's function Go is the usual diffusion 
propagator: 

Go{uJ,q) = [-tu + Dq^y\ (14) 

To calculate the disorder-averaged Green's function G in the presence of the random drifts 
one could develop a perturbation theory in terms of v, as we consider only the small disorder. 
However, the long-range character of the disorder correlations in Eq. (|TT]) (which is effectively 
present even in the model (c) as F{q) does not vanish when q —>■ 0) makes all the diagrams 
diverging as u; — >■ in the case d < 2. The situation is similar to that in the quantum- 
diffusion problem where the appropriate divergence is due to the interference effects. More 
convenient way to perform the ensemble averaging is to represent G as a functional integral 
over the conjugate complex fields ^(r) and (f{r): 

^^''"'^^ = — jT^WMvr^ — ' (^^^ 

where the effective action functional is given by 

S\(p, cp] = jdf'r [iuoTpLf + v^ (dcJp) (f — DdaTpda^] (16) 



The averaging is performed by the standard replica trick: the fields ^(r) and ^{r) and the 
functional integration in Eq. ([T5|) are iV-replicated and the independent averaging over the 
numerator and denominator in Eq. ([T5|) is justified in the replica limit A^ = (that should 
be taken in the final results). Taking into account the disorder both in the random drifts v 



and in the diffusion coefficients D which is defined by Eqs. (||) and (|TOD, one deduces the 
effective action 

S\<p, <p\ = {^0 + Sint + Scum} [^, f] (17) 



which should be substituted for that given by Eq.(|TBp into Eq.([T3|) where the functional 
integration should be performed over all components of the fields cp = (fi . . . (f^ and ip = 

ipi . . .^jv- Here 

So[ip, cp] = Jd'^rlpizLO + Dod^)cp (18) 

SU^,^] = y /ci'^rdV(9,^<p)^F„^(r-r')(9^W)^, (19) 
oo r ^ 

Scum [^, <fi]=^T. g^'^'5(«) ' ^(s) = d"^ n {dc.jpdc., <p) ^ (20) 



s=2 i=l 



Choosing the dimensionality of the fields cp and 'cp so that the action So is dimensionless, 
one sees that the random drift term, Eq. (|19[), is relevant in dimensionalities d < 2. The 
RG analysis of the functional (0), ([T9|) in the upper critical dimensionality d = 2 describes 
the anomalous long-time behaviour of the average transport coefficients. On the other 
hand, the naive counting of the scaling dimensions shows that the higher-order gradient 
operators, Eq. (pOD , are irrelevant in any dimensionality d > 0. Indeed, under rescaling 
L — > XL, g^*) — > A'-'^^-'^^'^g^''-' so that the naive dimension of g'-'^-* is always negative. Further 
on, we will show that, similar to the case of the quantum diffusion |T^, in the case of 
the potential model, Eq. ([TT|a) , the one- loop RG corrections overturn this conclusion and 
make the scaling dimensions of the high-order gradient operators positive which will lead to 
nontrivial asymptotic properties [§ of the distribution P{D). 



III. LONG-TIME BEHAVIOUR OF THE TRANSPORT COEFFICIENTS 



We begin with renormalizing the functional (p!8D , (p!9|) to find the long-time behaviour of 
the transport coefficients. The standard RG procedure is performed by expanding exp{iSint) 
in a power series, integrating with the weight exp {iSo\(pQ, c^q]) over the "fast" components 
of the fieldsand exponentiating the results of the integration. Here (p{r) (and ^(r)) is 
decomposed into the sum of "slow", ^(r), and "fast", V'o(^) components where 

^(r) = Y: <Pi^y''-\ ¥'o(r) = E '^(q)e^'^^ 

q<Mo A(jo<g<go 

< A < 1 is the scaling parameter, go is the ultraviolet cutoff, e.g. the inverse lattice 
constant in the lattice realization of the model. The actual (loop) expansion is carried out 
in powers of the dimensionless (at d = 2) disorder parameter 



7 
47rL)2 



(21) 



which takes a part of the effective couphng constant of the model. (Here we omit the 
subscript o iii 7 and D, and will keep it only to denote the bare, unrenormalized values of 
appropriate quantities). In the weak-disorder approximation considered here go -C 1. 

In the two-loop approximation, one comes in the replica limit (A^ = 0) to the following 
RG equations for the two parameters of the functional ([T8|), (|19D (as usual, the frequency 
uj is not renormalized due to the conservation of the total probability): 



dXnD/dhiX-' = ag - 2(1 



a 



2\„2 



(22) 



2\„2 



d\ii-i/d\n\-^ = -(1 - a)g - 2(1 - a^)g 
which lead to the following RG equation for the coupling constant 

dlng/dln\-^ = -{1 + a)g + 2(1 - a'^)g^ . 



(23) 



Here, for the models defined in Eq. (11 



a-c 



a 



-1, 
1, 
0. 



[a) 
(b) 
(c) 



(24) 



For the models of the solenoidal and mixed disorder, Eq. (^^,c), one faces the "zero-charge" 
situation: the coupling constant that is initially small (go ^ 1 in the weak-disorder case) 
is further decreased under the RG transformations, Eq. (^31). Therefore, the RG solutions 
of these model are asymptotically exact. In the potential-disorder model, Eq. ([TT|a), the 
coupling constant is not renormalized at all. This is obvious within the two-loop accuracy of 



Eq. (23). Moreover, this statement has been proved to the all orders of the loop expansion 
g] which, naturally, does not rule out a possibility of non-perturbative contributions to the 
coupling constant). The absence of renormalization of the coupling constant does not make 
the model to be trivial, as both the diffusion coefficient, D, and the strength of disorder, 7, 
do change nontrivially with the change of the scale A, according to Eq. (^2]). Solving this 
equation and setting the logarithmic RG variable to ln(DgQ/u;)^/^, or equivalently, in the 
time representation, to ln(t/r)^/^ (go is the ultraviolet cutoff which is of order the inverse 
lattice constant in the lattice realization of the model, and r ~ l/Dq^), one finds the 
long-time behaviour of the diffusion coefficient as 



D{t) 



d_ 

di 



r^{t)) oc < 



Do{t/T)~^^'^, potential disorder (a) 

-Do In^' ^(t/r), solenoidal disorder (b) 

Dq[1 + A/ Init/r)], mixed disorder (c) 



(25) 
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So, indeed, the potential and solenoidal models show sub- and super-diffusion behaviour, 
respectively, while the long-time asymptotics of D in the mixed-disorder model is practically 
unaffected. 

In a similar way, one finds the long-time dependence of the mobility, /i, by substituting 
the drift field in the presence of an external electric field E, V = ;uE + v(r) for v into the 
FP equation (^ and the effective functional ( [T8| ) and solving the appropriate RG equations. 



It yields 

;Uo(t/'r)^S/2^ potential disorder (a) 

/io, solenoidal disorder (6) (26) 

/io/ln(t/r), mixed disorder (c) 



/i(t) oc 



Note that the Einstein relation between the renormalized diffusion coefficient and the mobil- 
ity holds only for the potential model. Its violation for the other two models can be ascribed 
to the renormalization of the effective temperature []^. In the potential case temperature 
remains unrenormalized. Applying the results (^), ( PB| ) to the hopping model above and 
substituting g as in Eqs. (^) (P), one finds that the mobility (and thus conductivity) acquires 
a characteristic temperature dependence exp(— T~^). 

In general, all the three models show the scale-dependence of the averaged transport 
coefficients on the disorder which is totally different from that in the quantum-diffusion 
problem. In the latter, for d = 2 the effective disorder parameter, the inverse dimensionless 
conductance g~^, increases with the scale ITH] so that the system is driven towards the 



strong disorder (localization), while in the former the effective disorder parameter, g, either 
decreases with the scale (solenoidal or mixed disorder) or remains unchanged (potential 
disorder). This makes the classical-diffusion models considered to have asymptotically exact 
solutions but, unfortunately, prevents even qualitative considerations of the strong-disorder 
limit within the perturbative RG approach which proved to be so fruitful in the case of the 
quantum diffusion (see for reviews Ref. ||I2| ). 



Although the disorder-dependence of the average quantities in the models considered has 
nothing to do with that in the quantum diffusion problem, their mesoscopic properties show 
quite a striking resemblance to those of the quantum model. We demonstrate this, limiting 
further considerations to the potential-disorder case which is most interesting of the three. 



IV. MESOSCOPICS FLUCTUATIONS IN THE POTENTIAL-DISORDER 

MODEL 

In considering the "mesoscopic properties, we should answer for the two important ques- 
tions. First, whether the model considered shows the absence of self- averaging which is 

9 



characteristic for the quantum-diffusion problem at T = 0, so that the fractional value of 
the fluctuation does not vanish in the thermodynamic limit, L — > oo. We will show that this 
happens due to the long-range character of the correlations (|ll]a). Then, what is a natural 
"mesoscopic" scale beyond which the self-averaging is restored? In the quantum-diffusion 
problem, such a scale is a coherence length [0,0. There is no direct analog to that in the 
problem of non-coherent transport but the "mesoscopic" scale arises naturally within the 
lattice hopping model considered above. 



A. Mesoscopic scale 

In real lattices the correlation function (^ becomes non-singular (oc qaQp) for q ^ tq 
where tq is a screening radius. The continuum model with the correlation function (|ll]a) 
could then describe only the random walks at a distance not exceeding tq. In general, tq 
takes the part of the infrared cutoff similar to that of the phase-breaking length in the 
weak-localization theory and defines the mesoscopic scale for the classical-diffusion problem 
considered. The screening radius can be "mesoscopically" large for systems where the carrier 
density is much lower than the density of the charged impurities which are responsible 
for the asymmetric disorder in the hopping probability that leads to the presence of the 
random potential drifts in Eq. (Q). Among numerous examples are diffusion of a charged 



particle injected into a disordered dielectric or inversion layer |T^, or that in weakly doped 
semiconductors with a very high or a very low degree of compensation, K. In the latter 
case, the screening radius measured in the length of elementary hops is proprotional to a 
large parameter |]T^, (1 — K)~'^^^ for strong compensation (1 — K <ti 1), and K~^/'^ for weak 
compensation {K <^ 1). Note that in all the cases above disorder is not expected to be 
weak in contrast to our model. Nevertheless, by focusing on the influence of the asymmetric 
disorder which is dominant in the model, we hope to shed some light at the processes in 
real system. Our main aim, however, is to use the model as a toy one to demonstrate a 
possibility of deep similarity between coherent and non-coherent mesoscopics. 



B. The absence of self-averaging in conductivity 

We start with the continuum model with the unscreened correlations (]Tl|a). To make 
an analogy with the quantum- diffusion problem most transparent, we consider by way of 
example the fluctuations of the conductivity (Jafi- It can be defined in terms of the Green's 
function of the FP equation as 

(Ta(3 = O-Q l^ap ~ J2 J ^"('^) ^^^('^' ^' '^ ^) '^(^'' ^"' ^) ^^ '^^' ^^") ' (^7) 
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In the average, {(Jap) = p^Saf3 where p is the particle density. The fractional value of the 
conductivity fluctuations is determined by the correlation function 

/Ca/3;75 = (SaapSa^s) / CTq , (28) 

where Saap = <Jai3 — {(^ap)- In the lowest nonvanishing order of the perturbation theory, one 
obtains for cu — > 0: 

/C«/3;.. = [l) / (01 ^-^^^ [F.MF,M) + F^AcDF.M] ■ (29) 

The long-range character of the correlations (|ll]a) leads to a divergence in /C in the limit 
L — >■ oo that exactly compensates the L~^ factor attached to the integral: 

/do 
~~^ 1 (30) 



where Sap^^S = Sa/sS^yS + ^a^Sps + 

After the substitution 



(31) 



where g is the average dimensionless conductance, this expression exactly corresponds to 
that for the appropriate correlation function in the quantum-diffusion problem [l^. Not 



only the proportionality to the proper disorder parameter is the same but also both the 
tensor structure of this expression, and the dependence on a geometric shape of the sample 



which is hidden in the same diverging integral. Naturally, the scale beyond which Eq. (|30D is 
no longer valid is totally different for the two cases (the screening radius of charged impurities 
in the former, and the coherence length in the latter), as well as behaviour at larger scales. 
Nevertheless, the similarity is quite striking, and even more so in much more subtle properties 
of the fluctuations like the long-tail asymptotics of the distribution functions. 



C. Anomalous dimensionality of the high-gradient operators 

If the distribution of the fluctuations were Gaussian, only the average and the second 
moment would be relevant. So one should find the higher moments of the distribution 
to determine its shape. We show here following Ref. ||^ that the high moments of the 
distribution function becomes large as compared to the second one so that only the tails of 
distribution are distinctly non-Gaussian similar to those in the quantum-diffusion problem. 
The distribution of the diffusion coefficient is easiest for considerations as the couplings g*^^^ 
that define its higher moments are directly included into the effective action (|20|) . To find 
the scale-dependent moments, we perform the one-loop renormalization of these couplings. 
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We will not discuss this renormalization in detail (see Ref. 0) but outline the main steps. 
First, the index structure of the action ( pO[ ) is not conserved under the renormalization. A 
set of additional operators is generated having the structure 

d^-^'^dp^^d.^'^ds^"... (32) 

with all possible permutations of the vector indices a, /3, . . . and the replica indices a,b, . . ., 
each index being repeated twice, where the summation over repeated replica indices from 
1 to A^ and over repeated vector indices from 1 to d{= 2) is implied. These operators are 
"unphysical" in a sense that the distribution P{D) is defined by the renormalization of the 
initial operators (^) only. Then, only those operators are of interest which have the RG 
feedback to the initial ones. 

To classify all the additional operators one introduces matrix notations 

Q"' = d^-^'^d^^', V''' = d^ip''d^ip\ V''' = d^ip'^d^Tp' . (33) 



In these notations, the initial action (^OD is /(i'^r(TrQ)*. It is possible to show that the RG 
equations have a triangular structure: the operators containing only the matrices Q are not 
influenced on by those containing VV. Therefore, all the relevant operators may be written 
down as 

S,[Q] = ^(^) Jd'^r { [TtQY' . . . [Tr(Q)]''" . . .} . (34) 

Here the set of integers s = si . . . Sm ■ ■ ■ obeys the constraint 

J2 i^Sm = s, (35) 

m>l 

the initial operator ( pO|) corresponding to s = (0, . . . , 0, s). The renormalization of the s-th 
cumulants results from solving the RG equations for the whole set of g^^\ the bare values of 
all the additional charges being equal to zero. 

The couplings g^^^ may be represented as some ket-vector defined by the "occupation 
numbers" Sm 

g^^^ = \s) = \V'2'\..m'"^ ...). (36) 

The matrix set of the RG equations involving all the couplings (|36|) can be diagonalized 

exactly [^]. The eigenvectors are given by 

\p)=Y.9{s)xM^) (37) 

where the summation is performed over all the partitions {s} = si . . . Sm ■ ■ ■ of the integer 
s obeying the constraint (|35|) , g{s) = sl/nm^^^'^-^m' is the number of elements in the class 
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defined by the partition s, and Xp(s) are the characters of irreducible representation of the 
group of permutations characterized by the Young frame p having boxes of length pi . . . p^ ■ ■ ■ 
where \s)ummPm = s. The appropriate eigenvalues are given by 

a..(p) = filzl)+i: ''-<^--/'" + ^' (38) 

The maximum eigenvalue corresponds to the eigenvector characterized by the one-line Young 
frame with pi = s, pm = ioi m > 1 for which Xp(s) = 1 for all s, so that it is equal to s(s— 1), 
as in the case of the quantum diffusion. Note that it could be verified, without any reference 
to the representations of the permutation group, by mapping the renormalization group 
operator onto a certain one-dimensional model of bosons with nontrivial cubic interaction 
I^,|9[]. Thus, with the one- loop accuracy, the dimension of the operators coupled to the 



moments of the diffusion coefficient is given by 

as = -{s-l)d + gs{s-l), (39) 

so that for large enough s (s^g^^) the one-loop correction overtakes the negative naive 
dimensionality of the operators (^Ol). 

Note that there is a deep technical analogy with the quantum-diffusion problem. In the 
latter, the high-order moments of the conductance fluctuations are described with the help 



of the high-gradient operators in the nonlinear a model [14|. Their renormalization involves 



either the mapping onto the one- dimensional model of interacting bosons |14| or the analysis 



in terms of irreducible representations of the group of permutations [|19|, similar to the 
procedure outlined above, and leads to the anomalous dimension of the operators given by 
Eq. ( p9| ) after the same substitution (^) as for the variance of the conductance fluctuations 
(^). This results in nontrivial similarity between the properties of the fluctuations in the 
two systems. 

V. COMPARISON OF THE RESULTS FOR THE FLUCTUATIONS IN THE 
COHERENT AND NON-COHERENT DIFFUSION PROBLEMS 

There are two types of contributions into the fluctuations of the diffusion coefficient in 



the random- walks model considered, similar to the quantum diffusion problem [Q. The 
"normal" one is given only by the functional (pJO]) . It diverges in the infrared limit thus 
making the fractional fluctuation (((SD) )j / D"^ oc g^ to be independent of the size of the 
system, analogous to the UCF in metals. The appropriate contribution to the conductance 
fluctuation in the potential-disorder model is given by Eq. (|30D. The additional contribution 
to the diffusion cumulants ([T0|) is governed by the dimensions of the couplings g^''-' in the 
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high-gradient functional pO|). Keeping only the maximum eigenvalue, as in Eg. (pQ]) , one 
finds in the critical dimensionality d = 2 

{{{SDY)) oc (//L)2(«-i)-g^(^-i) , (40) 

where I is some microscopic length that could be of the order of the lattice spacing, etc. For 
s^g"^, these cumulants increase very fast with the system size L. It is quite similar to the 
"additional" contribution to the conductance cumulants (or density of states, or diffusion 
ones) in the quantum-diffusion problem which are proportional [Q to {I / L)"^ e^'^^'^~^^ with 
u = \n[go/g) where Qq is the value of the average dimensionless conductance g in the square 
of the size /^, and / in this case is the mean free path. To make the analogy more striking, 
one substitutes here the value of the parameter u in the weak-localization limit at d = 2, 
using g = Qq — ln{L/l) |]12|. It gives for the conductance cumulants in the quantum-diffusion 
problem 

(((<5<7)^))(x(//L)2(^-i)-^-^^(-i). (41) 



Therefore, the same substitution ( pT] ) relates not only the variance of the fluctuations in the 
coherent and non-coherent problems, but also the high moments of the fluctuations, Eqs. 
(0) and {^. 

The increase with L of the high-order moments leads to the lognormal asymptotic tails 
of the distribution functions |l^,p[] which are naturally identical after a proper deflnition of 
the parameters: 



f{6X) oc exp 



-InUsXc 
4?/, V 



___ a^ 

4u 



(42) 



Here 6X stands for either 6D or for g in the classical or quantum problems, respectively, 
a = L/l, and u is given above for the quantum problem, and equals gln(L//) for the classical 
one. 



VI. CONCLUSION 

We have shown that there exist a very deep and surprising similarity between the charac- 
teristics of the mesoscopic fluctuations in the conventional quantum-diffusion model and the 
model of the non-coherent ('classical') diffusion which is described in the continuum limit 
with the Fokker-Plank equation (H) with the quenched potential random drifts, Eq. (|ll]a). 
All such characteristics of one model may be obtained from those of another one with the 
help of the substitution (|3TD, i.e. by substituting a proper weak-disorder parameter. Such a 
parameter is deflned in a very different way for the two models. In the quantum-diffusion 
problem it equals to the inverse conductance while in the classical one it is proportional to 
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the inverse square of the diffusion coefficient (and, thus, to the inverse square of the con- 
ductance), Eq. (0). The similarity between the high moments of the fluctuations leads to 
the distribution of the diffusion coefficient in the classical model to be very similar to the 
conductance distribution in a weakly disordered metal. In both cases, the distributions turn 
out to be almost Gaussian in the weak-disorder limit but have slowly decreasing lognormal 
tails, and the part of the tails increase with increasing the disorder. The mathematical rea- 
son for the similarity is that the RG equations governing the cumulants of the distributions 
in both cases are classified according to the same irreducible representations of the group of 
permutations. However, the derivation of the RG equations for the high-gradient operators 
proved to be much easier in the classical-diffusion model than in the quantum one. 

This similarity occurs in spite of the fact that the average transport coefficients behave 
absolutely differently in the two models. Mathematically, this is due to a different behaviour 
of the coupling constants in the field-theoretical models describing the quantum and classical 
diffusion. The asymptotic freedom of the nonlinear a model that describes the quantum- 



diffusion problem pO[, i.e. the increase of the coupling constant (inversely proportional to 
the conductance) with increasing a scale, is believed to govern the Anderson transition. No 
transition occurs in the classical diffusion problem where in the case of the potential disorder 
a perturbative renormalization of the coupling constant proves to be absent in all orders |Q] 
thus leading to the sub-diffusion, Eq. (|25|a) . 

Then one can hope that it is possible to separate in some way description of the average 
values from that of the fluctuations. A very simple, if not too simple, assumption is a pos- 
sibility to use the one- loop (i.e. the lowest-order) RG results for the distribution, Eq. (^2|), 
by substituting more rigorous (or exact) results for the average quantities. Surprisingly, it 



allows to reproduce |IJ] exact one- dimensional results for lognormal distributions by substi- 
tuting the exact one-dimensional value of u into the formulae similar to (|4^) . It provides a 
basis for the conjecture that one can obtained a reasonable description of the distributions 
near the transition just by substituting e" ~ |f7 — gd'^ with a proper choice of the critical 
exponent u. 

Furthermore, there is a hope that studying the classical-diffusion problem described here 
gives a possibility to learn more about the mesoscopic properties of the quantum diffusion 
in disordered media. 
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